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Oh! Abstract 

We reconsider several current bounds on the variation of the fine-structure constant in 
models where all gauge and Yukawa couplings vary in an interdependent manner, as would 
^ ■ be expected in unified theories. In particular, we re-examine the bounds established by 
the Oklo reactor from the resonant neutron capture cross-section of 149 Sm. By imposing 
variations in Aqcd and the quark masses, as dictated by unified theories, the corresponding 
bound on the variation of the fine-structure constant can be improved by about 2 orders of 
magnitude in such theories. In addition, we consider possible bounds on variations due to 
their effect on long lived a- and /3-decay isotopes, particularly 147 Sm and 187 Re. We obtain 
a strong constraint on Act /a, comparable to that of Oklo but extending to a higher redshift 
corresponding to the age of the solar system, from the radioactive life-time of 187 Re derived 
from meteoritic studies. We also analyze the astrophysical consequences of perturbing the 
decay Q values on bound state /3-decays operating in the s-process. 



1 Introduction 



The nature of fundamental constants in physics is a long-standing problem. While certain 
constants can be thought of as merely unit conversions (c, K, k B , etc.), others such as 
gauge and Yukawa couplings can be thought of as dynamical variables. Indeed, such is the 
case in string theory, where the only fundamental parameter is dimensional, namely the 
string tension. The dimensionless gauge and Yukawa couplings are then set by ratios of 
the dilaton and moduli field vacuum expectation values to the string tension. Similarly the 
gravitational coupling (Planck mass) is scaled from the string tension by a modulus vev. 
Thus until these vevs are fixed, the fundamental coupling constants could vary in time. Of 
course, it is widely expected that non-perturbative effects will generate a potential for the 
moduli and fix their vev's (probable at some high energy scale), the mechanism and scale 
of this fixing are a subject of much debate. Thus in principle, one can consider variations 
in the fundamental couplings a logical possibility. 

Indeed, a considerable amount of interest in the possibility of time-varying constants 
has been generated by recent observations of quasar absorption systems. Observations of 
the energy level splitting between the 6*1/2 — ► P3/2 and S1/2 — > P1/2 transitions in several 
atomic states such as CIV, Mgll, and SilV, suggest a time variation in the fine structure 
constant by an amount Aa/a = (0.72 ± 0.18) x 10~ 5 ]l| over a redshift range of 0.5 - 3.5. 
In addition, there may be preliminary evidence for a variation A///// = (5.7 ± 3.8) x 10 -5 
in the ratio of the proton to electron masses fi = m p /m e for redshifts of ~ 2 - 3. 

Starting from the work of Bekenstein |||, there have been a number of attempts to 
formulate a dynamical model of a variable fine structure constant 0, [5]. These models 
typically consist of a massless scalar field which has a linear coupling to the F 2 term of 
the U(l) gauge field. The coupling of non-relativistic matter to the scalar field induces 
a cosmological change in the background value of this field which can be interpreted as 
a change in the effective fine structure constant. Independent of our prejudices (or lack 
thereof) regarding a fundamental theory, such models are difficult to construct in such a 
way as to remain consistent with other experimental constraints. For example, the presence 
of a massless scalar field in the theory leads to the existence of an additional attractive 
force which does not respect Einstein's weak universality principle. The extremely accurate 
checks of the latter || lead to a firm bound that confines possible changes of a to the range 
Aa/a < 10 -10 — 10~ 9 for < z < 5 H |7|, P in the context of the minimal Bekenstein model 
where a change in the scalar field is triggered by the baryon energy density. It was argued 
|| that a significant 0(1) coupling between the scalar field and the dark matter energy 
density is required in order to allow Aa/a ~ 10~ 5 and remain consistent with equivalence 
principle constraints ||. Thus it is natural to expect that in generalized Bekenstein models, 
not only the fine structure constant but all of the couplings and masses will depend on the 
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expectation value of a light scalar. 

In addition, there exist various sensitive experimental checks that coupling constants 
do not change (See e.g., ||). Among the most stringent of these is the bound on |Aa/a;| 
extracted from the analysis of isotopic abundances from the Oklo phenomenon [P|-||12|, a 
natural nuclear fission reactor that occurred about 1.8 billion years ago. While the Oklo 
bound |Aa/a| < 10~ 7 is considerably tighter than the "observed" variation, Oklo occurred 
at a time period corresponding to a redshift of about 0.14, and it is quite possible that 
while a varied at higher redshifts, it has not varied recently. That is, there is no reason 
for the variation to be constant in time. Big bang nucleosynthesis also provides limits on 



Aa/a [13, 14 1. Although these limits are weaker, they are valid over significantly longer 
timescales. 

The Bekenstein model and its modifications are introduced in ad hoc manner, and their 
relation to deeper motivated theoretical models is problematic. A major stumbling point 
on the path between the theory and phenomenology of a changing a is the masslessness 
of the modulus that mediates this change f[5]| . Indeed, to be relevant for the cosmological 
evolution now or in the recent past, the mass of this scalar has to be comparable or lighter 
than the Hubble parameter at z ~ — 5, whereas quantum corrections would tend to gen- 
erate a much larger mass. This is a generic problem for any interacting quintessence-like 
model that is similar to the cosmological constant problem. Since very little is actually 
understood about the latter, we do not think that this problem is a sufficient reason to dis- 
card phenomenological models of changing a. Disregarding the problem of masslessness of 
the modulus that renormalizes coupling constants, we proceed to analyze phenomenological 
constraints on a theory with a fixed (modulus-independent) high-energy scale Mp, unified 
values for all coupling constants at Mp, and a single modulus that changes the values of all 
coupling constants. Such a theory is motivated by a string model with dilaton-dependent 
coupling constants. One has to keep in mind, however, that the simplest string tree-level 
values for the couplings of dilaton to matter and rridiiaton = lead to a catastrophic non- 
universality in the gravitational exchange by this scalar, which violates the current bound 
by 10 orders of magnitude |5|, W\. One remedy to this problem may be a more complicated 
form of the dilaton- matter coupling with a universal extremum |TJ| . Another possibility is 
that a massless modulus contains a relatively small admixture of the string dilaton, so that 
all the couplings of this modulus to matter are suppressed to a level consistent with the 
equivalent principle |§. 

The possibility that significantly stronger constraints on the variation of the fine struc- 
ture constant can be obtained in the context of theories in which the change in a scalar 
field vev induces a change in the fine structure constant as well as the other gauge and 
Yukawa couplings was first explored in (see also, ||17|| ). There it was recognized that 
in any unified theory in which the gauge fields have a common origin, variations in the 
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fine structure constant will be accompanied by similar variations in the other gauge cou- 
plings. In other words, variations of the gauge coupling at the unified scale will induce 
variations in all of the gauge couplings at the low energy scale. Note that even in theories 
with non-universality at the string scale, there is almost always some relation between the 
couplings. 

It is easy to see that the running of the strong coupling constant has dramatic con- 



sequences for the low energy hadronic parameters, including the masses of nucleons [14 



Indeed the masses are determined by the QCD scale, A, which is related to the ultraviolet 
scale, Mjjvi by dimensional transmutation: 



a (M 2 ) = 9 ' Wv) - 4n (1 V 



where 63 is a usual renormalization group coefficient that depends on the number of massless 
degrees of freedom, running in the loop. Clearly, changes in g s will induce (exponentially) 
large changes in A: 

AA _ 2tt Aa s (M uv ) Aa s {M uv ) 
A 9a s (M uv ) a s (M uv ) a s (M uv ) ' 

where for illustrative purposes we took the beta function of QCD with three fermions. On 
the other hand, the electromagnetic coupling a never experiences significant running from 
Mjjy to A and thus AA/A 3> Aa/a. A more elaborate treatment of the renormalization 



group equations above Mz [M leads to the result that is in perfect agreement with [14 



AA Act 
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A a 

In addition, we expect that not only the gauge couplings will vary, but all Yukawa 



couplings are expected to vary as well. In |jT4| , the string motivated dependence was found 
to be 

— = ^ (1.4) 

h au 

where au is the gauge coupling at the unification scale and h is the Yukawa coupling at the 
same scale. However in theories in which the electroweak scale is derived by dimensional 
transmutation, changes in the Yukawa couplings (particularly the top Yukawa) leads to 
exponentially large changes in the Higgs vev. In such theories, the Higgs expectation value 
corresponds to the renormalization point and is given qualitatively by 

v ~ M P exp(-2-K c/a t ) (1.5) 

where c is a constant of order 1, and a t = h 2 /4ir. Thus small changes in h t will induce 
large changes in v. For c ~ h t ~ 1, 

Av Aau 

80 1.6 

v au 
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This dependence gets translated into a variation in all low energy particle masses. In short, 
once we allow a to vary, virtually all masses and couplings are expected to vary as well, 
typically much more strongly than the variation induced by the Coulomb interaction alone. 
Unfortunately, it is very hard to make a quantitative prediction for Av/v simply because 
we do not know exactly how the dimensional transmutation happens in the Higgs sector, 
and the answer will depend, for example, on such things as the dilaton dependence of 
the supersymmetry breaking parameters. This uncertainty is characterized in Eq. (|1 . 5| ) 
by the parameter c. For the purpose of the present discussion it is reasonable to assume 
that Av/v is comparable but not exactly equal to AA/A. That is, although they are both 
O(10 — 100) Aa/a, their difference |AA/A — Av/v\ is of the same order of magnitude which 
we will take as ~ 50 Aa/a. 

In [|14j] , these relations were exploited to derive a strong bound on variations of a during 
big bang nucleosynthesis. The standard limit |0| of Aa/a < 10~ 2 is improved by about 2 
orders of magnitude to Aa/a < 10~ 4 as recently confirmed in a numerical calculation [ HS] . 
Here, we will consider the effect of these relations on the existing Oklo bounds as well as 
derive new bounds relating to the long lived a- and /5-decaying isotopes, 147 Sm and 187 Re; 
we will also comment briefly on the influence of changing the fundamental couplings on 
s-process yields. 

Before proceeding, we note briefly that in the class of theories we are considering, we 
would predict that the proton-electron mass ratio is also affected. For example, we would 
expect that 

Aii AA Av 

~ — (1.7) 

[i A v 

From Eqs. (|1.3| ) and ( p..6|) , we estimate that A\i/\i ~ 3 x 10~ 4 , based on the reported 
claim of a variation in a jl| and is somewhat larger than that reported in 0. For related 



discussions see 20 



2 The Oklo Bound Revisited 

Approximately two billion years ago, a natural fission reactor was operating in the Oklo 
uranium mine in Gabon. Shlyakhter |9[] argued that a strong limit on the time variation of 
a was possible by examining the isotopic ratios of Sm in the Oklo reactor. This suggestion 
was confirmed by Damour and Dyson fllQ"| , who performed a detailed analysis of the isotopic 
ratios and the effect that varying a would have on the resonant neutron capture cross section 
of Sm. Their analysis provided a bound of 

< 10~ 7 (2.8) 



Aa 



a 
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The bound was derived primarily by calculating the shift in the resonance energy, E r ~ 
0.0973 eV, of the Sm neutron capture cross section which is induced by a variation in 
the Coulomb contribution. While a full analytical understanding of the energy levels of 
heavy nuclei is not available, it is nevertheless possible to obtain an estimate of the size 
of an energy shift if the fundamental parameters of the theory are varied. In particular, it 
is possible to identify the over-riding scale dependence of the terms which determine the 
binding energy of the heavy nuclei. We perform this estimate in the context of the Fermi 
gas model. We will argue that the Oklo data can provide a sufficiently strong bound on 
the variation of m q and A, or more precisely, on the variation of m q /A. 

We begin by considering the semi-empirical formula for the binding energy B(A, Z) of 



a spherical nucleus with mass number A and atomic number Z [21 



B(A, Z) =C V A-C S A 2/3 



a c Z 2 A- l ^ + 5(A,Z)a paiv A- 



+ C d Z 2 A- L + A shell (A,Z), (2.9) 



where the first three terms on right-hand side represent the volume, surface, and Coulomb 
contributions, respectively, and the last three terms represent corrections due to pairing, 
surface diffusiveness, and shell structure. The corrections due to asymmetry between neu- 



trons and protons are included in Cy and Cs- The coefficients in Eq. ( |2.9| ) are given 
by 
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Numerical values for the pertinent quantities are ay 



k = 1.79, a c = 0.717 MeV, a = 0.546 fm, r = 1.205 fm, a pair 

( 1 for even-A eve 

8(A,Z) = 



15.68 MeV, a s = 
11 MeV, e 



1 for even-A even-Z, 

for odd-A, 

1 for even-A odd-Z. 



(2.10) 
(2.11) 
(2.12) 
(2.13) 

18.56 MeV, 
0.5, and 

(2.14) 



The shell structure coefficients are discussed below. 

The Coulomb contribution has a simple interpretation as the electromagnetic energy 
stored in a uniformly charged sphere of total charge Ze and radius tqA^ 1 ^ 3 . The volume 
and surface contributions can be rewritten as: 



C v A-C s A 2 / 3 = -((T) + (V)), 



(2.15) 
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where T and V represent the kinetic and potential energy, respectively, of the nucleons. 
Based on the Fermi gas model and considerations of nucleon-nucleon interaction potential 
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<T> 
(V) 



5 (A- 2Z\ 
1 + 9(— ) 



Vv + Vsym 



A - 2Z\ 
A J 
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where 



3 kp 1 

~ Tom k Mri 



(2.16) 
(2.17) 

(2.18) 



is the zeroth-order contribution to the total kinetic energy and the terms with the coefficient 
fs represent surface correction. For a Fermi momentum kp = 1.36 fin -1 , To = 23.01 MeV, 



and the other quantities can be obtained as Vy 
fs = 1-184. 

Now consider the reaction 



-38.69 MeV, V s 



sym 



119 



Sm + n 



150 



Sm + 7. 



15.29 MeV, and 



(2.19) 



The Q- value of the reaction is 



Q = 5(150,62) - 5(149,62) 



-0.987Tb - 0.851VV - 0.245V sym + 1.615a c + 8.16 x 10~ 2 a pair - 0.172C, 



+A she u(150, 62) - A shell (149, 62), 



(2.20) 



where the numerical values for all the terms except for the last two on the right-hand 
side are -22.71, 32.93, -3.75, 1.16, 0.90, -0.21 MeV, respectively. Comparing the Q- 
value calculated above with the experimental value of 7.99 MeV gives an estimate of 
A s heii(150, 62) — A s heii(149, 62) = 0.33 MeV. Clearly, changes in T and Vy produce the 
largest effects on the Q-value. Thus not only will our limit be strengthened (relative to 
the purely electromagnetic bound fTO] ) due to the strong interaction, but also due to the 
enhanced sensitivity of the binding energy relative to the Coulomb term alone. 

Since we are lacking a complete analytical theory for the origin of the terms entering into 
Eq. (|2.9| ), it will be sufficient to concentrate on the dominant kinetic and potential terms. 
Even in this simplified formalism, it is not possible to account for the exact scaling of the 
dimensionful terms T , Vy and V sym . However, we can identify a certain degree of required 
scaling, and it is quite clear that an exact cancellation of such a scaling is extremely unlikely. 
Furthermore, the resonant cross-section proceeds through an excited state of 150 Sm, which 
happens to lie very close to the Q value given in Eq. (|2.20j ). Thus the quantity of interest 
is, 

E r = Q-E* = 0.0973eV (2.21) 
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where E* is the energy of the excited state of 150 Sm. Unfortunately, we expect that vari- 
ations in the fundamental couplings also lead to changes in E*. The previous bound on 
a was based on the presumption that taking into account the variation of E* with a only 
strengthens the bound, so that a conservative bound on Aa can be traced directly to AQ. 
Here we will have to rely on the probability that it is also highly unlikely that both Q and 
E* depend on all of the fundamental parameters in exactly the same way. We return to 
this point below. 

Before we derive our bound on possible variations of the gauge couplings it will be useful 
to first use Eq. ( |2.20| ) to derive the bound on a along the lines of [[H]]. If we ignore all of 
the unification arguments given in the previous section, then the only clearly identifiable 
piece with the electromagnetic coupling in Eq. ( |2.20| ) is the Coulomb term. Damour and 
Dyson (1996) |T(J derived the bound 

- 0.12eV < A{Q - E*) = AQ < 0.09eV (2.22) 

where 

AQ = Q° kl ° - Q now (2.23) 

For our purposes it will be sufficient to consider the limit \AQ\ < 0.1 eV. The Coulomb 
term in Eq. ( |2.20| ) is simply 

l.6a c ^ 1.16MeV ( —) (2.24) 
where «o is the present value of the fine structure constant. Thus 

AQ = 1.16MeV (-^j (2.25) 

We can therefore immediately derive the limit 

Aa 

a 

in good agreement with Damour and Dyson. 

We next attempt to use the same procedure to derive the limit on the gauge coupling 
when the unification argument of the previous section is included. We first note that if 
we simply associate all dimensionful quantities as originating from Aqcd, no significant 
limit is possible. In such a naive approach, one would argue that the masses of the light 
quarks can be neglected, and all hadronic parameters such as and the strength of the 
nucleon-nucleon potential scale linearly with A. It is clear, though, that in this limit there 
is no sizable effect on the position of the resonance, since 

AE r ~ A(Q - E*) ~ ^~ E ^ AA (2.27) 



< 1(T 7 (2.26) 
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Figure 1: The m q dependence of hadronic parameters, (a) illustrates the insertion of m q qq 
into the nucleon line that determines k n . (b) is the chiral loop that gives the m q dependence 
of the nucleon-nucleon interaction. 

and hence AE r /E r ~ A A/ A. Since the constraint on variation in E r is only 0(1) (cf. Eq. 
(|2.22|) ), one is simply left with A A/ A < 0(1). This point has been repeatedly stressed in 



the literature. For the most recent discussion, see, e.g. Ref. [[^ 



There are two generic problems that prevent a rigorous analysis of AE r as a function of 
m q /A. The first problem lies in fact on the interface of the perturbative QCD description 
and the description in terms of hadrons. In short, we do not know the exact dependence of 
hadronic masses and coupling constants on A and light quark masses. The second problem 
concerns modeling nuclear forces in terms of the hadronic parameters. 

Generically, the mass of a hadron i can be parameterized as 

mi~ const AM + Ki -± + 0{m%)\ , (2.28) 

where Kj reflects the dependence on the light quark mass. Clearly, masses of the members 
of the lightest pseudoscalar octet have a significant dependence on m q , and n n ~ 1/2. On 
the other hand, heavy hadrons remain massive in the chiral limit which suggests that for 
nucleons k n <C 1. In practice, Kjy (Fig la) probably varies from 0.1 to 0.2 due in part to a 
rather large value of the matrix element of m s ss over a nucleon state, (N\m s ss\N) ~ 100±50 
MeV |23[ , which is known to admittedly poor accuracy (the combined matrix element over 



u and d quarks is 45 MeV). Other particles which are important for nuclear dynamics 
such as vector resonances are expected to have kn < < K n . We note that although a 
complete analysis of Awij/Am, is certainly lacking (due to the difficulty of the problem), 
some insight can be gained from QCD sum rule analyses. 

Recent progress in understanding the chiral dynamics of few nucleon systems is not 
directly applicable to large nuclei such as Sm (See e.g., [^4] and references therein). Instead, 
one has to resort to a model description of nuclear forces such as contact interactions 
generated by t-channel exchanges of vector and scalar resonances. In this approach, the 
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relevant Lagrangian can be symbolically written as 



C = £ N{id^ - m N )N + —(NN)(NN) 



(2.29) 



where the sum runs over all relevant isospin and Lorentz structures. The scale M remains 
finite in the chiral limit, and thus scales approximately linearly with A. It is important to 
note that single pion exchange, which is the most important source of the nuclear force in 
the deuteron, plays little if any role in nuclei with large A. Nevertheless, one can provide a 
conservative estimate of the m q dependence of a nucleon-nucleon potential by considering 
chiral loop corrections to the contact nucleon-nucleon interaction, Fig. lb. This way one 
obtains 
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A 2 



1 



(N\m q qq\N) 
f* . 



j In 









+ 0{m 



2 ) 



(2.30) 



and f n scales linearly with A. In principle, the strength of overall numerical (A and m q 
independent) coefficient a can be fit to nuclear data but for the purpose of present discussion, 
we simply take it to be 0(1). Unfortunately, there is no known reliable way of calculating 
the coefficients a® and in front of the zeroth and first order terms of the expansion in m q . 
On the other hand, it is very unlikely that ao and a\ would conspire to form exactly the same 
combination of A and m q as in m^r, ( [2. 28 ). Fortunately, the presence of a chiral logarithm 
in (|2.30|) in the 0(m 2 ) term and the known absence of such a term in ensures that 
the nucleon mass, rriN, and the interaction strength, 1/M 2 , depend differently on (m q /A). 
Thus, to be very conservative, we consider only the variation of the chiral logarithm in 
(f2~29|) and use 100 - 200 MeV for the value of (N\m q qq\N)^, noting that the linear term 
(of order m q ) is likely to produce a more stringent bound: 



AC ~ (1 - 4) x lQ-Hm 2 (NN)(NN) 



A{m q /A) 
m q /A 



+ 



(2.31) 



From AC we obtain the Hamiltonian AH using the mean field approximation: 

A{m q /A) 



AH = (1 - 4) x 10 _1 fm 2 n 



(2.32) 



where R is a nuclear radius and n an average nuclear density inside the nucleus, n ~ 0.17 
fm -3 . The sum runs over all nucleons inside a given nucleus. As discussed above, the 
change in the position of the resonance comes as the result of the change in the position of 
the ground state of 149 Sm and the excited state of 150 Sm. In this naive approach, we have 
neglected all specific origins of the resonant energy found in Q and E* . If we then associate 
AE r with the difference in AH between the relevant states of 149 Sm and 150 Sm, we obtain 



\AE r \ (1 - 4) x lO-Wno 



E r 



E r 



A{m q /A) 
m q /A 



(2.5 - 10) x 10 7 



Av AA 

v A 



(2.33) 
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Note, that our result Q2.33|) is about an order of magnitude weaker than the result 
obtained in [2^]. Applying \AE r /E r \ < 1 [TIJ, we arrive at the following bound: 



Av 

v 



AA 



A 



< (1 -4) x 10" 



(2.34) 



In the theoretical framework discussed in this paper, the result ( 2.34 ) allows us to 
improve the limits on the variation of the coupling constant by over an order of magnitude 
compared to Ref. [10|. Indeed, combining eqs. Q2.34 ) and the estimate |AA/A — Av/v\ ~ 
50Aa/a, one gets 

< (2 - 8) x 10~ 10 (2.35) 



Act 



a 



We remind the reader that the range quoted in the limit above is due to the uncertainty 
in the strange quark contribution to the nucleon mass and corresponds to using 100 - 200 
MeV for (N\m q qq\N). 

It is also possible (though it bears its share of uncertainty) to use the expression for Q 
( 2.20|) to further strengthen the bound. One should note that the resonance energy, E r ^ is 
tiny due to a cancellation between Q and E*. Both are individually « 8 MeV. However, 
the value of Q is also determined by a set of cancellations between the terms in Eq. ( |2.20| ) 
which are of order 20-30 MeV. By once again relying on the improbability that all terms 
will depend on the fundamental constant in the same way, one can in principle use only the 
largest term (or terms) in Q, which are the potential and kinetic terms. 

If we again parameterize our ignorance of exact scalings with |AA/A — Av/v\ 
expect that 



we can 



AE r A(Q — E*) AV 



E r 



E r 



E r 



V_ 

E r 



Av AA 



A 



< 1 



which leads to 



Aa 



a 



< (1 - 5) x 10 



-10 



(2.36) 



(2.37) 



for the range 0.1 < k < 0.5. In Eq. (|2.37|) , we have assumed V ~ 30 MeV. This is the most 
optimistic bound that one can expect from the Oklo data. Recalling, that the Oklo event 
occurred some 2 Gyr in the past, we would obtain the limit a/at < 2.5 x 10~ 19 yr _1 . 



3 Constraints on the variation of fundamental cou- 
plings from long lived a- and /3-decay nuclei 

Bounds on the variation of the fundamental couplings can also be obtained from our knowl- 
edge of the lifetimes of certain long-lived nuclei. In particular, it is possible to use relatively 
precise meteoritic data to constrain nuclear decay rates back to the time the solar system 
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was formed (about 4.6 Gyr ago). Thus, for the current standard inputs of Q m ~ 0.35, 
f^A — 0.65 and ho = 0.71, we can derive a constraint on possible variations at a redshift 
z ~ 0.45 which borders the range over which such variations are claimed to be observed (z 
= 0.5 - 3.5). Note that nuclei with relatively short lifetimes, especially in the case of bound 
beta-decay (see section 4), could be of astrophysical interest in the context of s-process 
nucleosynthesis [p5|. 



3.1 Nuclear physics 

Here, we will concentrate only on the possible influence of variations in the Q-value on a- 
and /3-decay lifetimes. As fission processes are less well understood and less sensitive than 
the other decay processes, we do not consider those isotopes for which fission is the dominant 
decay mode. On the other hand, a- and /3-decay processes are better understood, and our 
basic constraints can be seen as due to a change in the phase space of emitted particles 
in the case of /^-decays, while for a-decays, the constraints are related to Coulomb barrier 
penetrability. 

In principle, slight changes in a coupling strength can stabilize (destabilize) certain 
isotopes. Clearly, the maximum effect is expected to occur for nuclei with small \Q\. A 
pioneering study on the effect of variations of fundamental constants on radioactive decay 
lifetimes was performed by Dyson |p6 |. The isotopes which are most sensitive to changes 



in the Q value are typically those with the lowest value of Qp, which is the Q- value corre- 
sponding to /3-decay. In addition, AQ's arising from variations of fundamental constants 
are expected to scale with AB, which is the change in the binding energy B of the par- 
ent nucleus. Accordingly, as a first step, we consider isotopes with the smallest values of 
\Qp\/ B as calculated from nuclear mass tables WA , regardless of their actual decay modes. 
For a-decay, small values of Q a (of order the smallest Qp, i.e. a few tens of keV) are not 
interesting because the Coulomb barrier is so high that the decay probability is vanishingly 
small. Hence, we considered isotopes with the smallest Q a values which have dominant 
a-decay mode (Q a ~ 1 MeV). 

The possible influence of the variation of constants on the binding energy and on the 
Q values has been discussed in section 2. Here, we will consider AQ/B as a parameter 
(assuming again that AQ scales with AB, i.e. AQ cx AB and AB/B oc Aa/a). Isotopes 
with small Q for /5 + -decay are not considered as the electron capture channel is already 
open and will dominate the decay. Tables 1 and 2 present respectively the a priori most 
interesting nuclei concerning /3~-decay and electron capture (EC). These isotopes have 
\Q\/B less than the typical value of 10~ 4 and are selected from the NUBASE files of 
nuclear data p7[ . For each isotope, the Q value, lifetime, binding energy B and the Q/B 
ratio are displayed; also shown is the difference in spin and parity between the parent and 
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Table 1: Properties of selected nuclei with low \Qp/B\ 



A 


Element 


Or- 

p 

fkeV) 


half life 

(Vear) 


decav 


B 

(MeV) 


GIB 


AI 


An 


106 


Ru 


40.0 


1.02 


0~ 
H 


9.07 x 10 2 


4.41 x 10" 5 


1 


+ 


107 


Pd 


33.0 


6.5 x 10 s 


0~ 

y 


9.16 x 10 2 


3.60 x 10" 5 


2 




123 


Sb 


-53.3 




stable 


1.04 x 10 3 


-5.11 x 10~ 5 


3 


+ 


150 


Nd 


-87.0 


2.1 x 10 19 


23~ 


1.24 x 10 3 


-7.03 x 10~ 5 


1 




151 


Sm 


76.7 


90. 


3~ 


1.24 x 10 3 


6.16 x 10" 5 







148 


Eu 


41.0 


0.15 




1.22 x 10 3 


3.36 x 10" 5 


1 




157 


Gd 


-60.1 




stable 


1.29 x 10 3 


-4.67 x 10~ 5 







160 


Gd 


-105.6 




stable 


1.31 x 10 3 


-8.07 x 10~ 5 


3 




163 


Dy 


-3.0 




stable 


1.33 x 10 3 


-2.26 x 10~ 6 


1 


+ 


171 


Tm 


96.4 


1.92 


P~ 


1.39 x 10 3 


6.96 x 10~ 5 







179 


Hf 


-110.9 




stable 


1.44 x 10 3 


-7.71 x 10~ 5 


1 


+ 


184 


Re 


31.5 


0.10 


P + 


1.24 x 10 3 


-7.03 x 10~ 5 


1 




187 


Re 


2.6 


4.35 x 10 10 


P~ 


1.49 x 10 3 


1.74 x 10~ 6 


2 





the daughter nuclei, which governs the degree of forbiddenness. Stable nuclei which could 
become unstable are also included. In Table 3, we show the a-decay isotopes with half-life 
longer than 10 6 years and 1.9 < Q a < 4.7 MeV. 

In the following, we concentrate on the most favorable cases and we study the variation 
of their half-life as a function of AQ/B in the limit of \AQ\/B < 10~ 4 . We see that the 
isotope with the smallest Qp- value is 187 Re (2.66 ± 0.02 keV). In the case of electron 
capture (Table 2) a few isotopes also have small \Q\ values, such as 187 0s and 163 Ho. These 
low Q isotopes could be involved in bound state /3-decay prevailing in stellar conditions 
(see section 4). 

The rate of beta-decay depends on the decay energy and on the leptons' angular mo- 
mentum £. The decays are classified as allowed (I = 0) and ^-forbidden (l ^0). Unique 
transitions occur when only one multipolarity is permitted by spins and parities. We will 
consider here only the allowed and first forbidden transitions and neglect those that would 
give an even longer half-life. The half-live, Ti/ 2 , is related to the nuclear matrix element, 
Mfi, through 

27T 3 ln2 

(m e ) 5 G 2 F \M fi \ 

where Gp = l/2v 2 is the weak coupling constant and / contains all the energy dependence. 



f T ^ = 7Z. ^2u, ,„ (3-38) 
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Table 2: Properties of selected nuclei with low Qec/B 



A 


Element 


Qe.c. 
(keV) 


half life 

(year) 


decay 


B 

(MeV) 


Q/B 


A T 

AI 


A 

Att 


107 


Ag 


-33.0 




stable 


9.15 x 10 2 


-3.60 x 10 5 


2 




123 


Te 


r 

53.3 


> 6.0 x 10 i4 


EC+ 


1.04 x ltr 


5.12 x 10 


3 


+ 


136 


Cs 


80.0 


3.61 x 10~ 2 


(3- 


1.14 x 10 3 


7.01 x 10~ 5 


5 


+ 


150 


Pm 


87.0 


3.06 x 10- 4 


P~ 


1.24 x 10 3 


7.04 x 10" 5 


1 




151 


Eu 


-76.7 




stable 


1.24 x 10 3 


-6.16 x 10~ 5 







157 


Tb 


60.1 


71. 


EC + 


1.29 x 10 3 


4.67 x 10~ 5 







160 


Tb 


105.6 


0.2 


(3- 


1.31 x 10 3 


8.07 x 10- 5 


3 




163 


Ho 


3.0 


4573. 


EC + 


1.33 x 10 3 


2.26 x 10~ 6 


1 


+ 


171 


Yb 


-96.4 




stable 


1.38 x 10 3 


-6.96 x 10" 5 







176 


Lu 


106.2 


3.78 x 10 10 


P~ 


1.42 x 10 3 


7.49 x 10~ 5 


7 




179 


Ta 


110.9 


1.82 


EC + 


1.44 x 10 3 


7.71 x 10~ 5 


1 


+ 


178 


W 


90.0 


5.92 x 10~ 2 


EC + 


1.43 x 10 3 


6.30 x 10" 5 


1 


+ 


184 


Os 


-31.5 




stable 


1.47 x 10 3 


-2.14 x lO" 5 


3 




187 


Os 


-2.6 




stable 


1.49 x 10 3 


-1.74 x 10^ 


2 




193 


Pt 


56.6 


50. 


EC + 


1.53 x 10 3 


3.70 x 10~ 5 


1 




194 


Hg 


40.0 


440. 


EC + 


1.54 x 10 3 


2.61 x 10~ 5 


1 




202 


Pb 


49.0 


5.25 x 10 4 


EC+ 


1.59 x 10 3 


3.08 x 10~ 5 


2 




205 


Pb 


51.1 


1.53 x 10 7 


EC + 


1.61 x 10 3 


3.17 x 10~ 5 


2 
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Table 3: Properties of selected nuclei with small Q c 



z 


A 


Element 


(MeV) 


half life 

(year) 


B 

(MeV) 


Q/B 


60 


144 


Nd 


1.905 


2.29 x 10 15 


1.20 x 10 3 


1.59 x 10~ 3 


62 


146 


Sm 


2.528 


1.03 x 10 8 


1.21 x 10 3 


2.09 x 10~ 3 


62 


147 


Sm 


2.310 


1.06 x 10 11 


1.22 x 10 3 


1.90 x 10~ 3 


62 


148 


Sm 


1.986 


7. x 10 15 


1.23 x 10 3 


1.62 x 10~ 3 


64 


150 


Gd 


2.809 


1.79 x 10 6 


1.24 x 10 3 


2.27 x 10~ 3 


64 


152 


Gd 


2.204 


1.08 x 10 14 


1.25 x 10 3 


1.76 x 10~ 3 


66 


154 


Dy 


2.947 


3. x 10 6 


1.26 x 10 3 


2.34 x 10~ 3 


72 


174 


Hf 


2.495 


2.0 x 10 15 


1.40 x 10 3 


1.78 x 10~ 3 


76 


186 


Os 


2.822 


2.0 x 10 15 


1.48 x 10 3 


1.90 x 10~ 3 


78 


190 


Pt 


3.249 


6.5 x 10 11 


1.51 x 10 3 


2.15 x 10~ 3 


90 


232 


Th 


4.082 


1.41 x 10 1U 


1.77 x 10 3 


2.31 x 10~ 3 


92 


235 


U 


4.678 


7.04 x 10 s 


1.78 x 10 3 


2.62 x 10~ 3 


92 


238 


U 


4.270 


4.47 x 10 9 


1.80 x 10 3 


2.37 x 10~ 3 



For allowed transitions, / is given by |28 | 

f(Z, R, e ) = J p(Z, R, e e ) (e - e e ) 2 e e \[e\ - 1 de e 



(3.39) 



where Z is the atomic number of the daughter nucleus, e (eo) is the (maximum) electron 
energy in units of m e , and 



p(Z,R,e e ) = 2{2p e R) 2 ^ 1 \l + s)e^ 



T(s + irj) 



r(2s + l) 



(3.40) 



is the relative electron density at the nuclear surface (radius R), with 77 = aZE e /p e and 
s 2 = 1 — [Za) 2 ■ We will use / to calculate the energy dependence of the half-life for the 
allowed and first forbidden transitions. For unique first forbidden transitions (AI = 2, 
An = — ), we will use the approximation /« ~ (e^ — 1)/ (see p9]| ). The calculated half-life 
is displayed in Figure |2] as a function of AQ/B for the isotopes of interest. The individual 
panels of Figure || correspond to \ AQ\/B < 10~ n for n = 3, 4, 5, and 8. 

When \AQ\/B is relatively large as in Fig. §a, we see that several isotopes may change 
from stable to unstable or vice versa. For smaller variations as seen in Fig. [|b, we see 
that the half-life is slightly altered for a few isotopes, but the most spectacular effects 
arise for 187 Re and 163 Dy, the latter becoming unstable. When we restrict variations to 
\AQ\/B < 10~ 5 , we see that only the 187 Re half-life is significantly altered. After an 
examination of Tables 1, 2, and 3 and the dependence of the half-life on AQ, we conclude 
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that Re is the most promising isotope for studying variations of fundamental constants. 
We will discuss 187 Re in more detail after considering the case of long-lived a-decays. 



3.2 Limits due to Long lived a-decays 

Dyson [^BJ defined the sensitivity of the decay constant of the nucleus to the change of the 
electromagnetic coupling constant as 

s =w (3 - 41) 

which is a function of the decay energy Q. The a-decay rate depends on a through the 
probability of Coulomb barrier penetration. Approximating this rate by 



A oc exp 



AiraZ 



^Q a /{2m p ) 

where Z is the atomic number of the daughter nucleus, we can write 



(3.42) 



AiraZ ( 1 , 

I - -u , (3.43) 



where 



QJ{2m p ) V 2 



Q a da 



The decay energy Q a is given by 

Q a = B(A, Z) + B a — B(A + 4,Z + 2), (3.45) 

where B a is the nuclear binding energy of the a-particle. Varying only the Coulomb terms 
in Eq. ( |2.9| ) and neglecting the contributions from B Q , we obtain s = 774, 890, 575, 659, 
571, 466, and 549 for 147 Sm, 152 Gd, 154 Dy, 190 Pt, 232 Th, 235 U, and 238 U, respectively. 

The variation Aa/a is related to AA/A as 

^ = i^. (3.46) 
a s A 

Since the values of s for the nuclei listed above are similar, the most stringent constraint 
on Aa/a is given by the nucleus with the smallest known AA/A. While it is tempting to 
use 238 U for this purpose, the uncertainty of 6.7 x 10 -4 given in Table 4, corresponds to 
a laboratory measurement and is not directly applicable to a constraint at high redshift. 
Indeed, 238 U is extremely well measured, and is used to calibrate the ages of the meteorites. 
Instead, we consider 147 Sm as an example, and assume that AA/A is less than the fractional 
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AQ/B(xl0 5 ) 

AQ/B (xlO 8 ) 



Figure 2: Evolution of half-life of the f3 -isotopes as a function of the AQ/B parameter. 
In a) the scale is taken to be AQ/B < 10~ 3 , in b) AQ/B < 1(T 4 ; in c) AQ/B < KT 5 , in 
d) AQ/B < 1(T 8 . 



16 



Table 4: Sensitive and insensitive a- and (3 -isotopes 



Nucleus 


Decay 


Q (MeV) 


Half-life (yr) 


AA/A lab (%) 


AA/A mct (%) 


238 TJ 


a 


4.27 


4.468 x 10 9 


0.067 




235 |j 


a 


4.678 


7.038 x 10 8 


0.071 




232 Th 


a 


4.082 


1.405 x 10 10 


0.4 


> 1 


147 Sm 


a 


2.31 


1.06 x 10 11 


1.9 


0.75 


87 Rb 


P~ 


0.283 


4.75 x 10 10 


0.8 


0.8 


4U K 


P~ 


1.311 


1.265 x 10 9 


0.6 


1 



meteoritic uncertainty of 7.5 x 10 3 in the half- life of 147 Sm |3(J given in Table 4 (see also 

Aa 



3TJ). This gives 



< 10" 5 . 



a 



(3.47) 



The change in A due to more general variations of the fundamental constants can be 
written as 

Aa 1 (Am p AQ C 



AA 
X 



AiraZ 



^Q a /{2m l 



a 2 \ m r 



Qa 



(3.48) 



We again take Sm as an example, for which 



Q a = 3.36T + 3.40VV - 0.070V; ym + 39.9a c - 0.975C d 
+A sh (143, 60) - A sh (147, 62) + B a , 



(3.49) 



where the numerical values for all the terms except for the last three on the right-hand 
side are 77.31, —131.51, —1.08, 28.63, —1.18 MeV, respectively. The experimental value 
of B a is 28.30 MeV. Comparing the Q a calculated above with the experimental value of 
2.31 MeV gives an estimate of A sh (143, 60) - A sh (147, 62) = -1.84 MeV. We again see that 
changes in T and Vy produce the largest effects on Q a . Note however, that because the 
decay-rate depends on the ratio of Q to m p , if we neglect quark mass contributions to Q 
and m p and retain only the scaling due to A, we see that the contributions from Am p /m p 
and AQ/Q in Eq. ( |3.48j ) cancel and no improvement in the limit is possible. Including the 
quark contribution through the coefficient, k, we expect that the scaling of T and Vy with 
m q /A will not be exactly the same as that of m p . So Eq. (|3.48|) gives 



AX 



2naZn T(V) A(m q /A) 



v /Q a /(2m p ) Qa K/A) 



(3 - 20) x 10 2 



Av 

v 



AA 



A 



(3.50) 



using A(m q /A)/(m q /A) = \Av/v — AA/A|. The range is due to the difference in T and V 
(80 - 130 MeV) compared to Q ~ 2 MeV and the range in k = 0.1 — 0.5 as discussed above. 
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The variation in Eq. fl3.50|) corresponds to 

Aa 

a 

for AA/A < 7.5 x 10~ 3 . As shown below, this constraint is less stringent than that derived 
from 187 Re /3-decay. 



< (0.8-5) x 1(T 7 



(3-51) 



3.3 Limits due to Long lived /3-decays 

In section 3.1, we have seen that 187 Re is the most sensitive indicator of a possible variation 
of a as first argued by Peebles and Dicke [|32| and Dyson [[26]. The upper limit obtained from 
the analysis of the Re/Os ratio in iron meteorites obtained at that time was a/a = 5 x 10~ 15 
yr _1 |[26|| , and was less stringent than the Oklo limit [^Oj] . In spite of this, Re is of interest 



since the estimate of the effect of the variation of the coupling constants (particularly when 
we go beyond variations in a) based on the resonance energy (in the Oklo case) is more 
complicated than that based on the Q value (in the Re case). As we saw in section 2, 
the constraints based on the Sm resonant energy required some knowledge of the role of 
quark masses in the nucleus and a limit based on variations of Aq C d alone could not be 
obtained. In the case of /3-decay, a constraint can be obtained in a more direct way. Above 
all, the Re analysis is independent of the Oklo analysis and uses different physics. Finally, 
the dramatic improvement in the meteoritic analyses of the 187 Re/ 187 0s ratio mandates an 
update of the constraint on variations of the coupling constants. 

Rhenium occurs in relatively high concentration in iron rich meteorites. The 187 Re 
decay rate has been determined through the generation of high precision isochrons from 
material of known ages, particularly iron meteorites. Using the Re-Os ratios of IIIAB iron 
meteorites that are thought to have been formed in the early crystallization of asteroidal 
cores, Smoliar et al. |53| (see also ||34|| ) found a 187 Re half-life of 41.6 Gyr within 0.5% 
assuming that the age of the IIIA iron meteorites is 4.5578 Gyr ±0.4 Myr which is identical 
to the Pb-Pb age of angrite meteorites [ 35] . For a general discussion see refs. [|36" , 3j| . The 



results of Smoliar et al. [[33| are in good agreement with those of Shen et al. ||37|| , which 
adds confidence to the meteoritic value of the half-life (which is more precise than the direct 



measurement |$8| which carries a 3% uncertainty). The ages of iron meteorites determined 
by Rhenium dating are in excellent agreement with other chronometers such as U-Pb and 
Mn-Cr, which means that the Rhenium lifetime has not varied more than 0.5% over the 
age of iron meteorites (4.56 Gyr). This gives A/A < 1.1 x 10~ 12 yr _1 , to be compared with 



the limit of 10 10 yr 1 given by Dyson PS|| . The improvement in the data is considerable. 

The /3-decay of 187 Re is a unique first fordidden transition, for which the energy depen- 
dence of the decay rate can be approximated as |39|] 

A oc G 2 F Qlm 2 e . (3.52) 
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which is in agreement with and gives a good description of the variation of T\ji with 
Qp shown in Fig. ||d. The decay energy, Qp, is given by 



Qp = 5(187,76) -5(187,75) + (m n -m p -m e ) 

= 0.339T + 0.61iy sym - 26.4a c + 0.807C d + A sh (187, 76) - A sh (187, 75) 

+ (m n -m p -m e ), (3.53) 

where m n is the neutron mass and the numerical values for all terms except for the last 
three on the right-hand side are 7.80, 9.34, —18.93, and 0.98 MeV, respectively. The 
experimental value of (m n — m p — m e ) is 0.78 MeV. Comparing the Qp calculated above with 
the experimental value of 2.66 keV gives an estimate of A sh (187, 76) — A sh (187, 75) = 0.03 
MeV. Considering only the variation of the Coulomb term in Qp, we have 



A Qp \ Qp I \ at J V a J 



which gives 5a /a < 3 x 10 -7 for 5X/X < 0.5% over a period of 4.6 Gyr or a/a < 6 x 10 -17 
yr _1 . This is 100 times more stringent than the constraint in [^6| due to the improvement 
of the limit on A/A. 

The contributions to Qp from T and V , which scale with A, are comparable to that 
from the dominant Coulomb term, which scales as a A. As changes in A are O(30) times 
larger than that in a, we can estimate 

^=3^-2^3^^2x10^, (3.55) 
A Qp v Qp A A 

which gives 

Aa n 

< 8 x 10~ 9 3.56 

a 

for AA/A < 0.5% over a period of 4.6 Gyr or a/a < 2 x 10 -18 yr _1 . 

Note that all of these limits based on Re decay hold only if the variation of A is of 
the same order as the accuracy of A. This hypothesis can be cross-checked by different 
chronometric pairs with different sensitivities to variations of a. Furthermore, we can check 
a posteori, that even though the meteoritic ages are determined in part by lab measurements 
of the 238 U lifetime, the limits above still hold. To see this, we note that the adopted 
uncertainty in the Re half-life is determined by the uncertainty in the slope of 187 0s/ 188 0s 
vs. 187 Re/ 188 0s. The uncertainty in the age of the meteorites is neglected. However, Re is 
far more sensitive to changes in a than is U (ie., the sensitivity factor for U is about 500, 
while for Re it is 2 xlO 4 ). It is relatively simple to check that a consistent limit requires 
changes in a which are sufficiently small so that the uncertainty in the meteoritic age can 
be neglected. 
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4 s-process nucleosynthesis 



A variation of the fundamental constants could have several other significant consequences 
on astrophysical processes particularly on nucleosynthesis. For example, variations in the 
gauge couplings would affect the position of the triple a-resonance necessary for the synthe- 
sis of 12 C. This has been examined recently by Oberhummer et al. [40|. Here, we focus on 
the nucleosynthesis of heavier isotopes and in particular those generated by the s-process. 
We will keep our discussion qualitative since the nucleosynthesis of neutron rich isotopes is 
complex at both nuclear and stellar levels. Branching on the s-process path occurs every 
time the /3-decay lifetime of a given isotope is commensurate with the neutron capture life- 
time (see ||41|| , for a review). Among the nuclei listed in table 1, several species are involved 
in s-process branching (specially in the Sm — Eu — Gd, W — Re — Os and H g — Tl — Pb 
regions, on the basis of their low Q/B ratio P2L H3, H4], s-process nucleosynthesis has 



been studied in the context of the classical constant temperature scenario and more realis- 
tically in connection to thermal pulses in AGB stars. In the constant temperature scenario 
of the classical s-process, thermal excitation of low lying energy states takes place, strongly 
modifying stellar lifetimes, and the subtle effects of the variation of the coupling constants 
are masked. However, more realistic models of the stellar s-process invoke highly convective 
situations during thermal pulses followed by long episodes of quietness. The AGB recurrent 
thermal pulses give rise to rapid mixing of freshly synthesized material in cooler zones. In 
these periods, the effect of variations in a could show up, since in the interpulse regime 
the thermal population of excited levels is suppressed. Thus, the situation is involved and 
deserves a dedicated analyses on the basis of refined nuclear networks coupled to realistic 
stellar models. 

There is however an interesting case which benefits from the high temperatures in- 
volved, that of metastable 180 Ta, which is the rarest isotope in nature and is assumed to 
be predominantly of s-process origin |EB, |£7|| , (see however 0]). It is also the only isotope 



which is stable in the isomeric state. Would it remain metastable if the coupling constants 
were different in the past? This question also deserves investigation. Here, we will assume 
that it does remain metastable. The difficulty in producing this isotope in contrast to 
the facility of its destruction by thermally induced depopulation of the short lived 180 Ta 
ground state in stars is reflected in its rarity. The population of an excited state of 179 Hf 
at relatively high excitation energy (7/2— at 214 keV) is crucial to the synthesis of the 
metastable 180 Ta. The decay rate of excited 179 Hf, in turn, is essentially dominated at high 
temperature by the bound state /3-decay, thus a small alteration of its Q value (say by a 
few keV) would have strong consequences on the final Ta yield (at least in the classical s- 
process context). As the Qp value for decay from the 214 keV excited state of 179 Hf is ~ 100 
keV, the variation of AQp/Qp would be limited to less than a few percent. Comparing this 
with AQp/Qp = (1/3)AA/A < 2 x 10 -3 for 187 Re, we expect that the constraints derived 
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from considerations of Hf decay would be ~ 10 times weaker than those presented in 
section 3.3 for 187 Re. 

More generally, bound state /3-decay is expected to occur in highly ionized media in 
which the decay electron has a high probability to be captured in an empty atomic orbit. 
It is the time reversed process of orbital (bound) electron capture. It occurs whenever Q 
(< 0) is of the order of the binding energy of electrons in the innermost shells. 187 0s with 
Q = —2.6keV is particularly sensitive to variations of the fundamental couplings. Other 
interesting cases are 121 Te, 163 Dy, and 205 T1. Thus, in stellar conditions, new disintegration 
channels could open up and stable nuclei in the laboratory could become unstable. Conse- 
quently, a slight perturbation of Q values could have significant consequences on the results 
of the s-process and more importantly on the Re/Os datation Indeed, the 187 Re decay 
rate could be considerably enhanced in stellar interiors [fPJ by the bound beta-decay of 
highly ionized 187 Re. At typical s-process temperatures (3.10 8 K), the bound state /3-decay 
of 187 Re in to the 9.75 keV 187 0s level is energetically possible, provided the degree of ioniza- 
tion is high. The non unique forbidden transition may give an overwhelming contribution 
to the 187 Re decay rate, non unique transitions being, in general much faster than unique 
ones. This effect could be easily suppressed by a slight change of Q related to a change of 
couplings. 



5 Summary 

We have considered the class of unified theories in which gauge and Yukawa couplings are 
determined dynamically by the vacuum expectation value of a dilaton or modulus field. 
While such theories may allow the possibility that the fundamental coupling constants are 
variable (in time), they general do so in an interdependent way. That is, one expects on 
quite general grounds that a variation in the fine structure constant is accompanied by 
a variation in all gauge and Yukawa couplings. Even more importantly, as described in 
nj, such variations are also accompanied by variations in quantities such as Aqcd and 



the Higgs expectation value which are determined from the gauge and Yukawa couplings 
by transdimensional mutation. Limits on the variations in these dimensionful quantities 
impose severe bounds on variations in the fine-structure constant. In the context of dy- 
namical models where the change in fundamenntal parameters is goverened by a nearly 
massless modulus, these limits are complementary to the constraints imposed by checks of 
the equivalence principle. 

Within this context, we have re-examined the constraints which can be obtained from 
the natural nuclear reactor at Oklo. The previous bound [10| of Aa/a < 10~ 7 was found 



by limiting the variations in the Coulombic contribution of the Q value for the resonant 
neutron capture process. We found that by including variations in Q which are induced by 
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variations in A and the quark masses, m q , we can improve this bound by approximately two 
to three orders of magnitude. The improvement is due to 1) the sensitity of both AA/A 
and Av/v to Aa/a (a factor of 50), and 2) the sensitivity to the dominant terms kinetic 
and potential terms in Q rather than the Coulomb term (a factor of 30). However, we 
lose a factor of 2 - 5, due to the rather uncertain contributions of the quark masses to the 
nuclear potential in a heavy nucleus. Thus we obtain the limit Aa/a < (1 — 5) x 10~ 10 . 
It is clear that this bound, which is valid at the time period corresponding to a redshift 
z ~ 0.15, would require severe fine-tuning in any model which attempts to fit the recent 
quasar absorption data with Aa/a ~ 10~ 5 over redshifts z = 0.5 — 3.5. In particular, our 
results strengthen the need for the fine-tuning in Bekenstein-type models, as emphasized 
ini. 

We have also considered the bounds which can be derived from long-lived or barely 
stable isotopes. Small variations in a (A and v ) can lead to changes in the Q value in 
heavy nuclei. We showed that limits from a-decay nuclei such as 147 Sm can be as strong 
as Aa/a < 10 -5 from purely Coulombic variations, and Aa/a < (0.8 — 5) x 10~ 7 based on 
147 Sm life-time uncertainties as small as AA/A < 7.5 x 10~ 3 . 

We further showed that improvements in meteoritic abundance determinations have 
enabled one to derive substantially stronger bounds based on the /3-decay lifetimes of 187 Re. 
From the analysis of the Re/Os ratio in meteorites with ages of 4.56 Gyr (known to an 
accuracy of 0.01%), the half-life of 187 Re has been determined to an accuarcy of about 
0.5%. Purely Coulombic variations in a lead to the bound Aa/a < 3 x 10~ 7 . From 
the age of the meteorites, this limit is applicable at redshift z = 0.45. Thus not only is 
it competitive with the Oklo bound numerically, but it corresponds to a higher redshift 
further emphasizing the difficulty in achieving a value of Aa/a ~ 10~ 5 at higher redshift. 
We further stress that the physics involved in deriving this bound is independent to that 
used for the Oklo bound. When variations in A are included, the bound is improved to 
Aa/a < 8 x 10~ 9 . 
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